Rhythmic tapping difficulties in adults who stutter: A deficit in beat perception, motor execution, or sensorimotor integration?

Objectives The study aims to better understand the rhythmic abilities of people who stutter and to identify which processes potentially are impaired in this population: (1) beat perception and reproduction; (2) the execution of movements, in particular their initiation; (3) sensorimotor integration. Material and method Finger tapping behavior of 16 adults who stutter (PWS) was compared with that of 16 matching controls (PNS) in five rhythmic tasks of various complexity: three synchronization tasks ― a simple 1:1 isochronous pattern, a complex non-isochronous pattern, and a 4 tap:1 beat isochronous pattern ―, a reaction task to an aperiodic and unpredictable pattern, and a reproduction task of an isochronous pattern after passively listening. Results PWS were able to reproduce an isochronous pattern on their own, without external auditory stimuli, with similar accuracy as PNS, but with increased variability. This group difference in variability was observed immediately after passive listening, without prior motor engagement, and was not enhanced or reduced after several seconds of tapping. Although PWS showed increased tapping variability in the reproduction task as well as in synchronization tasks, this timing variability did not correlate significantly with the variability in reaction times or tapping force. Compared to PNS, PWS exhibited larger negative mean asynchronies, and increased synchronization variability in synchronization tasks. These group differences were not affected by beat hierarchy (i.e., “strong” vs. “weak” beats), pattern complexity (non-isochronous vs. isochronous) or presence versus absence of external auditory stimulus (1:1 vs. 1:4 isochronous pattern). Differences between PWS and PNS were not enhanced or reduced with sensorimotor learning, over the first taps of a synchronization task. Conclusion Our observations support the hypothesis of a deficit in neuronal oscillators coupling in production, but not in perception, of rhythmic patterns, and a larger delay in multi-modal feedback processing for PWS.

The authors report a single study that compares timing performance of adults who stutter compared to controls for three tapping tasks that vary in rhythmic complexity. Also considered is the potential role of music training in mediating any differences between groups. A general assumption of the study is that there is larger tapping variability for adults who stutter compared to controls, with the goal of their study to identify / pinpoint the nature of the observed differences, using the Wing & Kristofferson (1973) model as one approach to decomposing tapping variability into separate clock and motor components. The main conclusion is that worse tapping performance of adults who stutter compared to controls is due to both increased central clock variability and increased motor variability.

Strengths
The topic of investigating potential timing deficits in stuttering warrants investigation, as there has been an increasing number of studies that have proposed timing deficits in stuttering, with mixed results. Studies that clarify the nature of timing differences between adults who stutter and controls are certainly needed.

Weaknesses
Although the topic of the study is interesting and warrants investigation, the reported study has a number of significant weaknesses that in my view preclude publication. It is further difficult for me to see how a revision would adequately address these concerns without effectively yielding a new manuscript. Here are the main issues.
• First, the methodology for the study is poorly described and motivated. The descriptions of the tasks are very hard to unambiguously interpret and also are non-standard versions of the task(s) to assess timing performance that have been used with the W&K model.
The section on task description (L250-299) has been written again to improve its clarity.
The introduction now details further the cognitive subprocesses involved in paced and unpaced tapping, in order to better introduce the research questions and to justify the different rhythmic tasks used in this study. Figure 1 gives a visual summary of the different tasks and the new Table 1 recapitulates which parameter was considered in each task.
With respect to the latter, it is not clear why the authors have chosen to use very non-standard versions of synchronize-continue tapping tasks, which makes it difficult to assess the data in relation to previous work. For example, the standard synchronization-continuation task has individuals synchronize finger taps with an isochronous series of tones, which after a certain number of tones, cut out, and participants continue tapping at the same pace/rate until there is a cue to stop. The continuation task is treated as a separate task -and seems as far as I can tell to be separate with participants listening during the synchronization phase and tapping during the continuation phase. In general, this alternative version of the sync-continue task is not well motivated and in general I would recommend using the standard version so that the data are more directly comparable to previous data.
First of all, we modified the label of that task to avoid confusion and comparison with a standard synchronization-continuation task, as used in many previous studies. The condition is now referred to as ISO_REPRO, since it aimed at testing the ability to perceive and reproduce a periodic pattern on one's own, at a specific tempo. The choice of starting by a passive listening phase, instead of a synchronization phase like in the classical synchronization-continuation paradigm, aimed at distinguishing 1-the ability to build an accurate representation of the beat from passive listening only, without motor engagement, and 2 -the possible "consolidation", or improvement in accuracy of that representation after several seconds of actual tapping. Our initial hypothesis was that PWS might present a tapping accuracy comparable to that of PNS immediately after passive listening, but that a group difference might emerge only after several seconds of taps, due to a deficient 'motor consolidation' mechanism in PWS. The choice of this non-standard task is now better explained in the material & methods section (L 265-274). The underlying theoretical background and the research question, in relationship to neural oscillations, are now more clearly presented in introduction (L108-122+ L194-213). The results section has been reorganized address these questions more explicitly (see L436-492). In particular, additional analyses have been conducted to test the possible improvement in tapping accuracy and consistency after motor engagement, compared to after passive listening only (L451-457 and L489-492).
The version [of the synchronization task] described in the manuscript seems somewhat odd as there are two tones presented (one high and one low), with the high tone indicating that the person should tap and tones are organized in eight element groups. Participants were then asked to tap for a certain number of repetitions of the 8-element pattern, but based on the text this was for only the first seven of the eight low tones (??).
The label of that task has also been modified to "1:1_ISO_SYNC" and its description has also been improved (L250-264) The choice of a metered isochronous pattern (instead of a simple metronome) was motivated by -Controlling for the metrical organization that would naturally arise when listening to an isochronous sequence (1-3), but that could vary from one listener to another. -Comparing a simple synchronization task with a more complex one, in which an external stimulus is played every 4 taps only, which clearly induces a perception of quadruple meter. It was therefore necessary to induce the same metrical organization in the first simple task, in which an external stimulus was played on every beat. -Exploring whether PWS show different abilities with meter compared with PNS, which might indicate a more global and lower level deficit in beat perception. There are further a number of confusing elements of the task and procedure descriptions. For example, in the description of the task order, the synchronization to complex rhythm task is listed twice and the continuation task is not listed at all. This seems like it has to be an error.
The description of the tasks has been improved and clarified (L250-299), and accompanied by the recapitulative Figure1. The label of the tasks has been changed in order to avoid the confusion between the complexity related to non-isochrony (in the condition now labeled "NONISO_SYNC") and the complexity related to hearing an external auditory stimulus every 4 taps only (in the condition now labeled "1:4_ISO_SYNC").
• Second, I find the description of theory in the introduction to be overly loose/muddled. In this respect, in the consideration of the view that individuals who stutter have a temporal processing deficit, the authors equate the theoretical concept of a central clock and with the generation of an internal beat.

This connection is problematic because what is meant by a central clock in the literature is a time interval measuring mechanism that is akin to 'stop watch' that measures the duration of each to-be-timed interval independently, with some variance (i.e., in the literature this is viewed as an interval-based timing mechanism) whereas internal-beat generation emphasizes the concept of a rhythm and a beat-based (entrainment) mechanism of timing. It is thus odd to be making connections between (correlated) estimates of central clock variance
where each successive interval is estimated independently and synchronization measures where synchronization relies on stimulus-driven entrainment. This feels like comparing 'apples' with 'oranges'.
Thank you for your relevant comment. The introduction has been completely re-written and complemented, following your recommendations. It now details the different subprocesses involved in paced and unpaced tapping, better defines the notions of "beat" and "meter" perception, and the theory of neural oscillators coupling (L108-122). The term "central clock" is now avoided, and the correlation between CCV (Central clock variance) and PLV (synchronization consistency) is no longer considered in the new version of the article.
• Third, there are a number of questions about how the dependent variables are measured or the resulting data itself, which reduces confidence in interpreting the results and drawing clear conclusions. As one example, the formula for measuring asynchrony seems like it has an error. As far as I can tell, the authors are taking the asynchrony in msec and then multiplying by 360 and then dividing by 0.5. Dividing by 0.5 is multiplying by 2, so they are taking the raw asynchrony in msec and then multiplying by 720, which would not yield a value between -180 and +180 degrees as claimed. The text also describes "with a 360 modulo", which does not make any sense to me. Here, I believe what the authors mean is that they took the raw time point of each tap (modulo 500 ms -the inter-beat-interval of the metronome), which if then divided by 500 ms would give a value between 0 and 1; they then would need to rescale between -0.5 and 0.5 and multiply by 360 to get a value between -180 and +180 (or an equivalent procedure). Note here also that the phase values of -180 and +180 correspond to the same phase, which is relevant for the calculation of mean phase angle. If the authors are simply averaging the phases to calculate mean phase angle, then average -180 and +180 would give a value of zero, which is incorrect. In order to compute mean phase angle, it is necessary to use circular descriptive methods.
Our measure of tapping asynchrony, as a phase angle is similar to previous studies (Sares et al. 2019;Falk et al. 2015). The value "0.5" corresponds to the Inter-stimulus Onset Interval (IOI, in seconds). So, dividing the asynchrony (between the stimulus onset and the tapping instant, also in seconds) by the IOI and multiplying it by 360 (degrees) gives a phase angle (in degrees Bayesian circular mixed models were already used in the previous version of the paper (R package bpnreg). Following your advice, PA results are now also presented with circular plots in the revised version of the manuscript. Nevertheless, we agree that it was problematic that the W&K decomposition was only valid for half of the tapping trains in the PWS group. We therefore decided to remove the analysis based on the W&K model and variance decomposition in the revised version of the article.

• With respect to the application of the W&K model to decompose tapping variability into clock and motor components, it is not clear whether the authors linearly de-trended the data before applying the model (removing drift), which is standard for use of this model. The model assumes a stationary time series of produced intervals and the effect of consistent drift is to reduce the negative correlation of successive intervals (reducing the estimate of motor variance -and sometimes yielding negative estimates). Positive correlation of successive intervals is a violation of the model and
• With respect to music training, it is not clear why dance was included as part of the music training measure. Dance is not music training, but an individual with > 5 years of dance and consistent practice would be classified according to the authors' procedure as having a high level of music training?
We considered that both the practice of dance and music can influence synchronization abilities, since they consist of producing gestures in synchrony with external rhythms. It has been proven that amongst musicians, drummers and professional pianists show a particularly high synchronization accuracy and consistency, better than singers or nonmusicians (Krause, Pollok, and Schnitzler, 2010) (4). Other studies also showed greater synchronization abilities in skilled dancers, compared to non-dancers (Miura, Kudo, Ohtsuki, and Kanehisa, 2011) (5). In any case, the 2 dancers in our study also had a corresponding level of musical training. So, to avoid confusion or discussion, we removed any mention to dancing in the assessment of Musical Training (see section A of the supplementary material). Furthermore, Musical Training was no longer considered as an explicative factor in the revised version of the article. It was only used as an individual factor, like gender and age, to match participants between the PWS and PNS groups. • A final comment about dependent measures is that the authors indicated that they used a peakpicking algorithm in MATLAB to determine peaks, but provide very little additional detail about how this algorithm determines peaks, which can be tricky. Along these lines, did the authors conduct any low pass filtering of the data before picking peaks to remove noise? Details of the peak-picking procedure need to be spelled out so that the reader can better evaluate the method used.
More details on peak detection are now given on L328-334: "First, the force signal was low-pass filtered (Chebyshev filter, cutoff frequency of 100 Hz, using the function filtfilt in Matlab (R2018b) to extract its envelop, and normalized, based on its maximum value observed in each executed tapping task. For each tap, the first sharp peak of the force signal, corresponding to the tapping instant, was detected automatically (using the Matlab function "findpeaks", with a minimum interpeak distance of 200ms and a 20% threshold for peak height). These tapping instants were saved in PRAAT (73) annotation files, and were all manually verified and corrected." • Fourth, I have two general comments about the results that detract and significantly limit the contribution of the work. It is unclear (and problematic)  Following your advice, as well as those of the second reviewer, we no longer consider Musical Training as an explicative factor in the revised version of the article. It is now only used as an individual factor, like gender and age, to match participants between the PWS and PNS groups. The figures, the results section and the discussion have also been modified accordingly and do not mention that factor any longer.
• Along these lines, I find it odd that although one of the main conclusions of the study is about the relation/correlation between estimates of central clock and motor variance during continuation tapping with the synchronization measures, none of the graphs show these correlations. Rather, Figures 6, 7, and 8, focus on group comparisons of PWS and PNS for the three levels of music training -which given the very small sample sizes for each level of music training are not very meaningful and unrelated to the main question of interest.
As indicated above, we removed any mention to musical training. All the correlations tested in the article are now supported by graphs and presented in section B of the supplementary material.
• Finally, I do not find the level of precision in the writing be up to par for publication. The writing would need to be significantly improved and some sections rewritten to improve precision and clarity to a level that is publication quality.
Many sections of the manuscript have been rewritten. Furthermore, we had the manuscript proof-read by a native English speaker. We therefore hope that it now reaches the expected level of precision and clarity.

Reviewer #2:
The The section Material & Methods has been rewritten to improve its clarity. The tasks have been renamed and are summarized in the Figure 1. The extracted measures are now summarized in the Table 1.

There were numerous naming inconsistencies throughout, for example, line 222, it is unclear what is meant by "metronome and tapping instants". Is this all tasks?
It is now clarified (L330-333) that tapping instants were detected from the first sharp peak of the force signal, recorded with the force pressure sensor. It is now also clarified (L257-260) that external auditory stimuli, or triggers, were made of both a metronome click marking the pulse, and of audio beeps marking the rhythmic pattern to reproduce. The composition of these external auditory stimuli in each rhythmic condition is summarized in Figure 1.

In the discussion (ln 614): greater tapping variability during unpaced tapping is mentioned: is this referring to the synchronization continuation task?
The discussion has been fully re-written. Greater tapping variability in PWS has been found in our study as well as most of the previous studies, for both synchronization tasks and tasks in which participants tapped on their own, without any external auditory reference (continuation tasks, reproduction tasks like ISO_REPRO here, or self-paced tasks).
Continuation isn't a pure unpaced tapping measure as they had a cue to begin. I was also wondering why there was no pure unpaced tapping task, as people who stutter have been shown to be aided by an external cue (i.e., as they had in the CONT task at the beginning).
True. Following your comment, we now avoid the term "unpaced" throughout the article and talk, instead, of tapping without the help of external auditory reference.

Inter-response interval could be more clearly labelled as inter-tap interval to fit with previous literature. These inconsistencies and presentation really need to be improved otherwise it's very difficult to follow the results.
Done. IRI was changed to ITI.
• Musical training One of the main concerns I had while reading these analyses was with the statistics related to musical training. From Table 1, there are only 3 participants in the "moderate" music training group and 3 participants in the "high" music training group in each stuttering vs. non-stuttering group. All the analyses including musical training are therefore reflecting very few participants, with a big group difference compared to those with "low" music training (n = 10 in each group). These different group sizes are also covered up by bar graphs, and it's impossible to see the spread of data, and whether there are important outliers. Musical training is also confounded with stuttering severity, as for the "moderate" group, there was 1 very mild stutterer and 2 mild stutterers, and in the "high" group there was 1 very mild, Following your advice, as well as those of the first reviewer, we do not consider any longer Musical Training as an independent variable of our experimental design. It is now only used as an individual factor, like gender and age, to match participants between the PWS and PNS groups. The figures, the results section and the discussion have also been modified accordingly and do not mention that factor any longer.

• Analyses
Were there convergence issues in your linear models? Adding musical training and severity as categorical fixed factors (with three levels each) into your model seems like it would have lots of problems, considering e.g., there are only 3 participants with moderate or high music training, and within each group, for those with high training, 1 is very mild, 1 is mild, and 1 is severe. It doesn't seem like you have enough data to model these interactions, and I would assume that R will tell you this. The statistical analysis section (2.7) seems to suggest that you could combine all of these factors (musical training + severity) in one model, but I couldn't find this in the results themselves.
The different mixed models considered are now clearly presented at the beginning of each paragraph in the Results section. Musical Training is no longer considered as an independent variable. Correlations with stuttering severity (though the SSI score) are tested in the PWS group only, and independently from the effect of other factors, such as Group, Task or Beat strength.

• Figures Individual variation should be displayed in all graphs by including individual data points, and/or better representations of the spread of the data (e.g., box plots, but individual data points would be ideal). This would allow the reader to easily see the spread/variance of data, and also group size differences between bars.
Following your advice, the figures now represent the average value observed for each participant (N=16 in each group), instead of standard deviations with errorbars.

• Discussion
The numerous theories presented in the discussion also make for some tough reading, with no strong conclusions being made. It seems in the end that it's unclear what the results show and how they could be reflected in the different models. Perhaps a clearer summary or more integration across theories is necessary here. The final conclusion that "the dual premotor model and the sensory accumulation model" are compatible with most of the observations didn't come out easily from the discussion. Some reframing and streamlining seems necessary here.
The theoretical framework of the whole article has been improved and complemented. Instead of considering two hypotheses: that the reduced tapping accuracy and consistency of PWS may be related to a timing deficit vs. a motor deficit, we now consider different subprocesses that may be involved in sensorimotor synchronization and that may possibly be impaired in PWS (1-motor execution, 2-beat perception and reproduction, 3sensorimotor integration and learning). The whole article is now organized around the examination of these abilities. This also enables clearer conclusions on the cognitive levels that appear unimpaired, and those at which a significant difference is observed between PWS and PNS. Thus, the results from our study, point towards (1) a deficit in neural oscillators coupling in production, but not in perception, of rhythmic patterns in PWS, and (2) a larger delay in multi-modal feedback processing for PWS.
• Paragraph starting line. 658 -starts suggesting that there was support for a global deficit in motor skill. I therefore expected this paragraph to show this. However, the conclusion of the paragraph is that stuttering is NOT caused by differences in motor skill. Please make a topic sentence that is consistent with the evidence presented in the paragraph.
The discussion section has been fully rewritten. It now makes a clearer distinction between 1-motor difficulties in movement initiation, 2-other possible sources of increased motor variability, and 3-possible deficits in sensorimotor control and learning.
In that framework, our results support the idea that PWS present (1) a deficit in neural oscillators coupling in production, but not in perception, of rhythmic patterns and (2) a larger delay in multi-modal feedback processing • Minor Points: Abstract: Authors mention that there are three finger-tapping synchronization tasks, but then they list 5. Figure 1 also lists 5. It would be useful throughout to be more consistent with the labelling of each task and order of presentation, to make it easier for the reader to process.
There are indeed three synchronization tasks: -a simple synchronization task with a 1:1 quadruple metered isochronous pattern (1:1_ISO_SYNC) -a synchronization task with a quadruple metered isochronous pattern, where only the strong beats (every four taps) are marked by an external auditory stimulus (1:4_ISO_SYNC) -a synchronization task with a non-isochronous pattern (NONISO_SYNC) and two other tapping tasks, where participants do not synchronize with auditory stimuli: -a "reaction" task, in which participants follow as quickly as possible an unpredictable pattern (REACT) -a reproduction task, in which participants reproduce on their own a quadruple metered isochronous pattern, just after listening to it passively (ISO_REPRO) Pg. 4, lines. 66-69 -Can you explain the Wing and Kristofferson method, or rephrase the sentence so the reader isn't expecting an explanation? Is there a reason it can only be applied on unpaced tapping (ln. 75)?
Since the W&K decomposition could not be applied to a substantial part of our data, this analysis was removed from the article and no mention is made to it any longer.
Pg. 5, lines 106-107 -Couldn't central clock variance be related also to motor execution problems?
The whole idea of the W&K decomposition was to disentangle, from the total tapping variance, what can be attributed to Central Clock Variance (CCV) and what can come from Motor Implementation Variance (MIV). We decided to approach the data differently and remove the analysis from the article, as well as from the discussion of its results. The discussion has been fully rewritten.
109 -could tapping force just measure confidence?
Yes, this is possible, and could be used in future studies to evaluate indirectly the difficulty of a task, for example. In this study, it was measured seaking 3 objectives: 1-to test whether participants, and in particular PWS, were able to perceive and reproduce the quadruple meter of the proposed patterns 2-to test the sensory accumulation theory, i.e. whether the degree of NMA correlates with the tapping force 3-to have a descriptor of the amplitude variability of finger movements, and not only of their timing variability Data Cleaning: were any taps excluded from the analysis? E.g., while they were beginning the task? From section 2.6 line 215 it seems that all taps were included? Could this increase variability?
Following your comment, the revised version of the article now distinguishes the very first taps in a task (first 8-beat cycle) and the "stabilized" phase (2 nd and 3 rd 8-beat cycles). Instead of simply excluding the very first taps, we found it relevant, in the two conditions 1:1_ISO_SYNC and ISO_REPRO, to explore the variation in tapping accuracy and consistency over time, between the very first taps and the stabilized phase, and to discuss these variations in terms of motor engagement and sensorimotor learning.
Pg. 15, were CV, CCV, MIV, IRI, Finger RT all in the same model? Aren't there big collinearities between these measurements? And if there's only 55% of the PWS group with MIV and CCV calculations, it's missing a lot of data (The MIV and CCV estimations were considered in the analysis only in these cases, which represented 68% of the tapping trains (82% of the PNS group and 55% for the PWS group) and no single value could be calculated for 4 PWS participants.) We agree that this was problematic that the W&K decomposition could apply to only half of the tapping trains in the PWS group. We therefore decided to remove any mention of the W&K model and variance decomposition in the revised version of the article. The article still presents variations in RT, RT_Var = std(RT), mean(ITI), PE = abs(0.5 smean(ITI)), and CV = std(ITI)/mean(ITI). No significant correlation was observed between CV and RT_Var (see section B.2.2 of the supplementary material). No significant correlation was also observed between PE and RT (see section B.2.3 of the supplementary material). The variation of each parameter was explored independently, with different mixed models, as clearly indicated at the beginning of each paragraph in the Results section.
Line 318 -was there a reason not to use the Watson-Williams test here?
A Watson-Williams test unfortunately did not enable to consider the factor "Participant" underlying the repeated measures of a same individual. Applying a Watson-Williams test over the whole set of "aggregated" data would not be correct, since it would artificially inflate the cohort size (as if we had 32*N repetitions of participants, instead of 32 participants). On the other hand, applying a Watson-Williams test over the 32 average value of each participant would be correct, but would considerably reduce the statistical power of the analysis. Instead, using mixed models (and here Bayesian circular mixed models) is the most recommended method to deal with repeated (angular) data, since it enables to consider the whole dataset, and to account at the same time for the factor "Participant" underlying the intra-individual variation of the dependent variable.
Line 325 -why is a generalized linear model suddenly used here? What distribution was used?
A generalized linear model was used, when testing the correlation between two variables that were measured in different conditions, or when testing the correlation between a variable and the SSI score, since repeated measures could not be considered in that case.
Only average values measured for each participant in a condition could be compared, which therefore did not require to use mixed models and to consider a random factor on the participant.
In the revised version of the article, correlations (that are presented in section B of the supplementary material) are simply tested with a Pearson's correlation test for linear data, and with an angular-linear correlation test (with the toolbox Circstats) for PA.
Lines 331-335 -then we have coefficients bc and SAM -it's unclear what this adds to the analysis.
In the first version of the manuscript, the reporting of these 3 parameters followed the recommendations of Cremers & Klugkist (2018) for the use of the bpnreg package and the reporting of correlations, based on Bayesian mixed models. As indicated above, the correlation of PA with other linear variables is now simply tested with an angular-linear correlation test (with the toolbox Circstats) in the revised version of the article. Figure 2: Can you please show individual data in this graph (e.g., as small dots)? Please also mention how many participants are in each group.

Results
Done. The size of each group (N=16) is also given in the legend.
The phase angles would be better represented as a circular plot in my opinion, e.g., by using the library "circular" in R, or in Matlab using the CircStat toolbox.
PA variations are now represented with circular plots, following your recommendation. Figure 1: from Ln 365-365 it seems that there were both strong and weak beats in SYNCSimp -can you include this information in the figure? Was there an emphasis placed on these beats? Otherwise how are they considered as strong?
The description of that task, now labelled 1:1_ISO_SYNC, has been clarified L252-260, and clearly indicates that: "Since a metrical organization of beats (into groups of 2, 3 or 4) arises naturally and automatically when listening to an isochronous sequence of identical tones (1-3), we controlled for that perceptual grouping and induced the perception of quadruple meter, i.e. with a "strong" or accentuated beat sensed every four pulses, the other beats sensed as "weak" or unaccentuated). To achieve this, auditory stimuli were organized into 8-beat cycles, with a metronome click marking the pulse on each beat, and an additional audio beep (Pitch: 1100 Hz; 20 ms) played simultaneously on the first seven beats only (without variations in pitch, loudness, or duration) No emphasis was placed on the 1 st and 5 th beat of the stimuli, which we reckon, is disputable and could be improved in a next study. However, we verified that participants perceived that quadruple metrical organization, since they tapped with greater force on these 1 st and 5 th beats (see section 3.3.2 and Figure 7).
The sheer number of acronyms in the results makes it almost impossible to follow at times.
We improved the notation consistency throughout the article, and gave a summary of the task and variable names in the Figure 1 and the Table 1. We hope that it is now easier to follow.
Line 561: Please outline again what REAC means, or use consistent terminology so it's clear which task is which.
REACT stands for the "reaction" task. The description of the tasks has been improved and summarized in Figure 1.

Typos:
Pg. 5 line 91, observed should be observe Line 95 -"these evidences of" should be "this evidence for" Line 717: in this line of "though" Line 760, has two commas. Please fix up others throughout These typos have been fixed or corresponded to paragraphs that have been fully re-written.